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• Abstract 

O ! 

, Theoretically motivated smallness of the penguin amplitude in B —* ttit decays allows 

to calculate the value of the unitarity triangle angle a (^2) with good accuracy. The 
Q^- relatively large branching ratio of the decay into 7r°7r° is explained by the large value of 

, FSI phase difference between decay amplitudes with 1 = and 1 = 2. 

^ PACS: 12.15.Hh, 13.20.He 

^ 1 Introduction 

The exclusive decay amplitudes of hadrons are determined by dy- 
namics at large distances and can not be calculated with presently 



available QCD tools. 



O 

Fortunately it was found long ago that the experimental data on 
branching ratios and CP asymmetries of B — > tttt decays allow to 
determine the value of the unitarity triangle angle a with essentially 
no hadronic input using isospin invariance of strong interactions only 
|IJ. However, large experimental uncertainties in particular in the 
values of the direct CP asymmetries lead to poor accuracy in the 
value of a determined in this way. 

If the penguin amplitudes are negligible in charmless strangeless 
B decays we would determine the value of unitarity triangle angle a 
from CP asymmetry S+- extracted from B — ► 7r + 7r~ decay data with 
essentially no theoretical uncertainties. As it was found in paper [2j 
neglecting penguin amplitudes one gets the values of angle a from CP 
asymmetries in Bj decays to 7r + 7r~, p + p~ and 7r ± p =F consistent with 
the global fit of unitarity triangle. Since the penguin contributions to 
these decays are different j2| the fact that the numerical values of a 
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are close to each other testifies in favor of smallness of penguin am- 
plitudes. Small penguin corrections to these decay amplitudes were 
accounted for in j3j where the hadronic amplitudes were found from 
the quark amplitudes with the help of factorization. However, it is well 
known that the branching ratio of Bd(Bd) —> tt°tt° decay predicted by 
factorization appears to be more than 10 times smaller than the ex- 
perimental data. The way out of this contradiction could be large FSI 
phases in B — > tttt decays. The validity of this theoretical ingredient 
will be checked by the more accurate experimental data. 

Though the penguin contribution is relatively small compared to 
tree amplitudes and can be neglected in the first approximation in the 
decay probabilities and in the CPV parameters S it determines the 
CPV parameters C and should be accounted for in the analysis of the 
complete set of observables. 

The charmless strangeless B decays are described by b — > uud quark 
transition. The effective Hamiltonian responsible for this transition 
consists of two parts: the tree level weak amplitude (operators 0\ 
and O2 in standard notations) dressed by gluons and the gluon pen- 
guin amplitudes (operators O3 — Oq); the parametrically small elec- 
troweak penguins are omitted. The gluon penguins being very im- 
portant in AS = 1 strange particles nonleptonic weak decays are 
almost negligible in AB = 1, AS = transitions. The reason is 
twofold: firstly, Wilson coefficients are much smaller in case of B 
decays because infrared cutoff is at fi ~ instead of fi ~ ^qcd] 
secondly, the enhancement factor originated from the right-handed 
currents m^/m s (m u + m^) ~ 10 for strange particles decays is re- 
placed by m^/m^mu + rrid) ~ 1/3 for beauty hadrons. That is why 
after presenting the general phenomenological expressions for the am- 
plitudes we will start our analysis of B — > tttt decays in Section 2 by 
the sequestered Hamiltonian which does not contain penguin contri- 
butions . From the experimental data on Bd(Bd) — > tt + tt~, 7r°7r° and 
B u —* tt + tt° branching ratios we will extract the moduli of the ampli- 
tudes of the decays into tttt states with isospin zero Aq and two A2 and 
find the final state interaction (FSI) phase shift 5 = 62 — So between 
these two amplitudes. The value of the unitarity triangle angle a in 
this approximation is directly determined by CP asymmetry S+-. 

-"^Let us stress that while from the smallness of B — > tt°tt° decay width it would follow that penguins are 
small, the opposite statement is not correct: the relatively large width to neutral pions does not necessary mean 
that penguins are large. 
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While the absolute values of the amplitudes Aq and A2 are re- 
produced with good accuracy by the factorization formulas, the FSI 
phase shift appears to be unexpectedly large, 5 = —(53° ±7°). This is 
the reason why B — > 7r°7r° decay probability is significantly enhanced 
in comparison with the naive factorization approach, where one ne- 
glects 5. In Section 3 FSI phase differences in K — > tttt, D — > tttt 
and B — > Z)7r decays are considered. In all these cases the phases are 
large, which is attributed to the existence of / = resonances in tttt 
scattering in the cases of K — > tttt and (partly) in D — > tttt decays 
while large FSI phases in B — > Dtt decays are unexpected. In Sec- 
tion 4 we consider the theoretical estimates of 5 and show how FSI 
can enhance B width to neutral pions not enhancing that to neutral 
p mesons in accordance with experimentally observed suppression of 
B — > p°p° decay width. 

In Section 5 the penguin contributions are considered; the cor- 
rections to the numerical values of Aq and 5 due to gluon penguin 
amplitudes are determined, as well as the correction to the unitar- 
ity triangle angle a and the values of CP asymmetries C + _ and Coo- 
In Conclusions the pattern of the B — > tttt decay amplitudes emerg- 
ing from the experimental data is presented. Appendix contains the 
calculations of the decay amplitudes in factorization approximation. 

2 B — > tttt without penguins: decay amplitudes 
from branching ratios 

The quark Hamiltonian responsible for B — > tttt decays has the parts 
with AI =1/2 and AI = 3/2 which produce 7r-mesons in the states 
with 1 = and 1 = 2 correspondingly. QCD penguins having 
AI =1/2 contribute only to the 1 = amplitude. Taking into ac- 
count the corresponding Clebsch-Gordan coefficients and separating 
the penguin contribution (P) with the CKM phase different from that 
of Aq we obtain: 
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where are CKM matrix elements and the penguin amplitude with 
an intermediate c-quark multiplied by V u bV* d + V c bV* d + V t bV td = is 
subtracted from the penguin amplitudes with intermediate it-, c- and 
t-quarks (the so-called t-convention) 2 . To check if the factorization 
works in B — ► 7T7t decays it is convenient to introduce /+(0) - the 
value of the formfactor which enters the amplitude of semileptonic 
B d —> it I is decay at zero momentum transfer in Eqs. (l)-(3). 7 and 
(3 are the angles of the unitarity triangle; 6 2 and 5q are FSI phases 
of the tree amplitudes with 1 = 2 and 1 = (below we will use 
6 — 62 — 60), 6 P originates from the imaginary part of the penguin 
loop with c-quark propagating in it p] while 6q is long distance FSI 
phase of the penguin amplitude. 60 in general is different from do] in 
Section 4 we will argue that pp intermediate state generate large value 
of do while its contribution into 60 is smaller: (pseudo)scalar part of 
penguin operator do not produce p mesons. 

The charge conjugate amplitudes are obtained by the same formulas 
with substitution (3, 7 — » —j3, —7. 

The CP asymmetries are given by [6j: 
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1-|A 7 



1 + |A 7 
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where tttt is tt + tt~ 



Or 7T°7T°. 



2 We prefer t-convention since the penguin contribution to the amplitude Aq in it is suppressed as 
m c/( m b mTO wV TO b)j while in the c-convention one subtracts the penguin amplitude with intermediate i-quark 
thus making the penguin contribution to Aq comparable to the penguin term. 
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Now we have all the necessary formulas and neglecting the penguin 
contribution we are able to determine Aq, A 2: S and the value of the 

unitarity triangle angle a from the experimental data on , Bqo, 

B + q and which are presented in Table 1. By definition: 

B + EE l/2[Br(B d TT+tT) + Bl(B d -> 7T+7T 

5 00 = l/2[Br(£ d 7r°7r°) + Br(5, 7r°7r° 

5 +0 = Br(fl u - tt+tt ) = Br(B u - tTtt ) , 

the last equality holds as far as the electroweak penguins are neglected. 

To extract the product A 2 f + {0) from B + q we will use the value of 
\Vub\ obtained from the general fit of the Wolfenstein parameters of 
CKM matrix (CKM fitter, summer 2005): A = 0.825 ± 0.019 , A = 
0.226 ±0.001, p = 0.207 ±0.040, fj = 0.340 ± 0.023 : 

\V ub \ = (3.90 ± 0.10) -10- 3 . (4) 

From (PJ) and the experimental data on B + q from the last column 
of Table 1 we readily get: 

A 2 f+(0) = 0.35 ±0.02 . (5) 

In order to understand if the factorization works in B u — > 7r + 7r° 
decay we should determine the value of /+(0). We find it using the 
data onB^ wlu decay from jBj: 

/+(0) = 0.22 ±0.02 , (6) 

thus getting: 

A 2 = 1.60 ±0.20 , (7) 

which is not far from the result of factorization: 

8 2 
A( = — t=(ci ± c 2 ) = -y=( ai ± a 2 ) = 1.35 , (8) 

see Appendix. We come to the same conclusion as the authors of 
paper [9j: A 2 is estimated correctly by factorization. Neglecting the 
penguin contribution we are able to extract the values of Aq and FSI 
phases difference S from Eqs. (CEJ)-(J3I) and the experimental data for 
£> + _, Bqq and B + q from the last column of Table 1. In this way we 
obtain: 

A) = 1.53 ±0.23 , (9) 
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which should be compared with the result of factorization: 



4 = ^(5ci - c 2 ) = 1.54 , (10) 

see Appendix. In this way we come to the conclusion that factorization 
works well for the moduli of both decay amplitudes. 
For the phase difference 5 = 62 — So we get: 

V3 B + - - 2i3oo + l%B +0 

COSd = — ; ; t= = , (11) 

4 J?B +0 jB + - + Boo-&B + o 



T+ -+uy- + - 1 ~uu 3t+ - 

£ = ±(53° ±7°) , (12) 

where tq/t + = r(Bd)/r(B u ) = 0.92 is substituted. This is the place 
where the factorization which predicts the negligible FSI phases fails. 

In Section 4 we will present a model in which the pattern of B — > tttt 
amplitudes obtained above is realized. 

Let us turn to the bottom part of Table 1. Since we neglect pen- 
guins the experimental value of S + - is directly related to the unitarity 
triangle angle a: 

sin2a T = S+_ , (13) 
"babar = 99° ± 5° , ^10 = 111° ±6°, aL r a g e = 105° ± 4° , 

where index "T" stands for "tree" stressing that penguins are ne- 
glected (three other values of a are not compatible with the Standard 
Model). 



3 FSI phases in K —> tttt, D —> tttt and B —> Dtt 

The s-wave amplitudes of two pions production with 1 = and 1 = 2 
are generally different. In particular there are quark- antiquark reso- 
nances in s-channel with 1 = but not with 1 = 2. This can lead to 
large difference of phases in the channels with 1 = and 1 = 2. Let 
us remind what experimental data tell us about these phases at the 
pion center of mass energies E = tjik and E = tud- Since at E = tjik 
only elastic rescattering of pions is possible (the inelastic channels 
are closed since the energy is low) Watson theorem is applicable and 
strong interaction phases of matrix elements of K — > (2tt)i decays are 
equal to the phases of amplitudes describing tttt — > tttt scattering at 
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E = ttik- From the analysis of tttt — > tttt scattering data performed 
in (ID] at E W7T = ttlk we have: 

^ = 35° ±3° , (14) 

8$ = -7° ±0.2° , (15) 

8$ -8$ = 42° ±4° . (16) 

The large value of 8q is due to the specific behaviour of / = 0, J = 
7T7r-phase attributed to the /o(600) (or a) "resonance". 

The same value of the difference 8q — &f follows from the analysis 
of Ks — > 7T + 7r~ , Ks — > 7T°7T° and K + — > 7T + 7T° decay probabilities anal- 
ogous to one we perform for B — > tttt decays in Section 2 neglecting 
the penguin contributions. (In case of K — > 7nr decays the penguins 
are very important being responsible for the enhancement of I = 
amplitude. Since CKM phase of the penguin amplitudes is almost the 
same as that of the tree amplitude, the analysis performed in Section 
2 is applicable for kaon decays but the amplitude Aq should contain 
the penguin contribution as well.) 

What concerns the moduli of the kaon decay amplitudes with 1 = 
and / = 2, they are given with rather good accuracy (within 50% from 
the experimental data) by factorization jllj . 

In case of D — > tttt decays, the gluon penguin amplitudes are negli- 
gible in comparison with the tree ones (since the loop with s-quark is 
subtracted from the one with d-quark while the momentum transfer 
is of the order of m%), and the effective Hamiltonian responsible for 
these decays looks like: 

H D = —j= sin c [q u^ a {\ + 75) defy* (1 + 7 5 )c + 

+ C?<fy a (l+75)d«7a(l+75)c] , (17) 

where 9 C is Cabibbo angle, sin# c = 0.22, cf « 1.22, cf sa -0.42 jEJ. 

Calculating the matrix elements in the factorization approximation 
we obtain: 

M D ^ +7T - = ^ S mQ c f°(0)f % m 2 D (c° + ^) , (18) 

= %Bin0j£(O)/X>(y + c?) , (19) 
= ^sin^(0)-^m^(cf + eg) , (20) 
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(the analogous formulas for B decays are derived in Appendix) while 
the isotopic analysis gives: 

M D ^- = ^smej?(O)Uml(0=A° + -±=Aj>) , (21) 

M D ^o n o = ^w0jD(p)f vm %(^A? - -j=A°) , (22) 

M D± ^o = ^sm6J°(0)f n ml(^^A°) . (23) 

From the recent study of semileptonic .D-meson decays D — > nlu it 
was found [T3j : 

/^(0) = 0.62 ±0.04 . (24) 
Comparing (J2BD, (EH) and recent measurement (Hj: 

Br(L» ± t^tt ) = (1.25 ± 0.10) • 10~ 3 (25) 

we obtain: 

A? = 0.88 ±0.08 , 
which is not so different from the factorization result, ( 1201 ): 

(A?) f = ^=(c? + c?) = l.2 . (26) 

Comparing Eqs. JED- JED with Eqs. J23J-® we obtain in the 
factorization approximation: 

(A°) f = ^J(5cf - e?) = 1.8 , (27) 



while according to |14j from the experimental data it follows: 

A? = 2.1 ±0.2 . (28) 

We see that the factorization results are within 30% from the exper- 
imental values of the moduli of the decay amplitudes. However, fac- 
torization fails completely in describing the difference of FSI phases. 
The data onD^ ir + 7r~ , D —> 7r°7r° and D ± —> 7r ± 7r° branching ratios 
lead to [Hj: 

6g-6l? = 6 D = ±(86° ±4°) , (29) 

which is responsible for (or follows from) the relatively large D — ► 7r°7r° 
decay probability [Hj: 

Bt(D° -> 7r°7r°) cxp = (0.79 ± 0.08) • 10~ 3 , (30) 



of the order of that into tt + tt (E] : 

Br(L»° 7T + 7T-) cxp = (1.39 ± 0.07) • 10~ 3 . 

Using td q /td + = 410/1040 from [THJ we readily reproduce the phase 
difference given by with the help of (JTTJ) . 

In factorization approximation neglecting ^ we will get: 



Br(D° - ttV)/ = i 



3 C 1 ~r <^2 
3 V c l + c 2 



2 

X 



xBr( J D+ ttV 3 )^ « Br( J D° 7rV) exp . (31) 

r D+ 

The analogous phenomena we encountered in B — ► 7T7t decays. 

Let us note that the s-wave resonance with zero isospin /o(1710) 
alone cannot explain such a big phase; its contribution to D — > 
(7T7r)/ = o decay amplitude is proportional to: 

m D - m fo -i^ = 150 MeV - *70 MeV , |^| « 30° . (32) 

It is not easy to reconcile reasonable (20%-j-30%) accuracy of factor- 
ization in describing the moduli of the decay amplitudes into 7ttt states 
with a definite isospin and the large FSI phases difference since the 
latter signal of strong rescattering of pions at E w tub which should 
not only generate phases but also shift the moduli of the amplitudes. 
The resolution may be that the interactions are "semistrong" in both 
channels: one half of 86° comes from / = 0, another from 1 = 2 (just 
as in the case of B — > tttt decays, see Section 4). 

If we suppose that FSI phases scale with decaying meson mass as 
1/M we will get about 30° phases difference for B —> tttt decays from 
(29). 

tttt FSI phase shifts at E = rrix and E = ttid are not small. How- 
ever, in both cases we are in the regions where two pion resonances 
are situated, which is not the case for the high energy of the order of 
5-meson mass. 

Our last example is B — > Dtt decays, where the energy is high 
and we are definitely above the resonances domain, though the FSI 
phase shift is nevertheless large 3 . Dtt pair produced in 5-decays can 
have / = 1/2 or 3/2. From the measurement of the probabilities of 



3 We are grateful to A.E. Bondar who brought this case to our attention. 
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B -> D°7r-, B° -> and 5° Z}°7r° decays in paper jEj the 

FSI phases difference of these two amplitudes was determined: 

5 Dn = 30° ±7° . (33) 

Concluding this section we wish to note that the direct CP asym- 
metry observed in Bd(Bd) ~^ tt t K ± decays is incompatible with small 
FSI phase difference between 1=1/2 and 7 = 3/2 amplitudes. 

4 FSI phases: theoretical considerations 

There are many theoretical papers on the final state interaction (FSI) 
in the heavy-meson decays [Ej - [22| . For example in paper [TBj the 
final state interactions in B — > tttt decay are modelled as the soft 
rescattering of the certain intermediate two-body hadronic channels 
(7T7T, pp, D*D*, DD). The hadronic amplitudes, which enter the calcu- 
lation of the imaginary parts of the decay amplitudes were described 
by 7r, /9, D, Z}*-meson exchanges in the ^-channel. Rather large phases 
due to FSI have been obtained. While for 7r-exchange (with the pole 
close to the physical region) this procedure is reasonable, it exagger- 
ates the contributions of the vector exchanges (p, D*) : which for the 
elementary particle exchange with spin J = 1 gives the partial wave 
amplitude, which does not decrease with energy. In reality all the ex- 
changes should be reggeized and in the physical region of the processes 
corresponding intercepts c^(0) < 1 (for D* the most probably value 
of the intercept is negative). This will lead to the strong reduction of 
the corresponding amplitudes (see for example [19|). 

A number of papers Q19J - [22J) use Watson theorem in order to 
extract the phases of the decay amplitudes by multiplying the bare 
matrix elements by S 1 ^ 2 (where S is the S-matrix) . However in B- 
decays there are many coupled multiparticle channels. In this case 
such a procedure can be applied only in the basis of the eigenstates 
which diagonalize ^-matrix. But for the realistic strong interactions 
this is impossible at present. 

Another approximation is to use the Feynman diagrams approach 
taking only the low mass intermediate states X, Y into account. This 
approach coincides with the use of the unitarity condition only if the 
transitions 7T7t — ► XY are described by the real amplitudes. This 
is certainly not true for elastic 7T7r-scattering, where the amplitude 
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at large energies is predominantly imaginary. In this formalism the 
resulting decay matrix elements are: 

ML = (6 vX S vY + iT J x f^) , (34) 

where M X y are the decay matrix elements without FSI and T^y ^ is 
the J = partial wave amplitude of the process XY — > 7T7t 4 . At very 
high energies the amplitudes of tttt elastic scattering are imaginary and 
TxY-^im d° n0 ^ decrease with energy (mass of a heavy meson). Thus 
this contribution, according to ( IH41) does not change the phase of the 
matrix element, but only changes its modulus. The extra phases come 
from the real parts of the amplitudes, which in Regge model are due 
to the secondary exchanges (R = p, /, ...), which decrease with energy 
as l/s 1-Q! fl( ) « l/yfs. The contribution of the pion exchange in the 
t-channel, which is dominant in the process pp —> 7T7t 5 decreases even 
faster (as 1/s). However Br (B — > p + p~) is substantially larger than 
Br (B — ► tt + tt~) and pp intermediate state is important in ( 1341) . This is 
especially true for color suppressed B — > 7r°7r° decay, where the chain 
B — > p + p~ — >• tt°tt° is enhanced. On the contrary tttt contribution 
in B — ► decay is relatively suppressed. Using Regge analysis of 
tttt scattering j23j and 7r-exchange model for pp — > tttt transitions, 
we obtain phases due to final state interactions for tttt final state 
62 ~ —12° and Sq « 18°. Thus the phase difference ^ « —30° is 
generated by intermediate pp and 7T7T states 6 . 

Sign of (5 is negative, just as in the case of K — > tttt decays. In this 
way in the numerical estimates we will use negative value of 5 from 

flU: 

6 w -(50 ± 7)° , J = 30° , £ 2 = -20° . (35) 

As far as Bi(B — > p + p~) is much larger than Br(5 — ► tt + tt~) 
because of enhancement in tree amplitudes but not in penguin am- 
plitudes (contribution of penguins in B — ► pp amplitude is small) we 
should expect 6q to be substantially smaller than 5q. 

Note that in this model there is little change in moduli of ampli- 
tudes in comparison with factorization predictions. 

For B — > pp decays the same model gives w — 5° for 7 = 2 and 
+5° for 7 = amplitudes, resulting in a small phase difference ~ 10°, 

4 We use the standard normalization with T J = s 2 ~ 1 . 

5 In B-decays transverse polarizations of p-mesons are small that is why a2 and oj exchanges in pp — > tttt 
amplitudes are suppressed. 

6 An accuracy of this number is about 15°. 
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consistent with experiment. 

Thus the lowest mass hadronic intermediate states may produce the 
phases which are consistent with the data on B — > 7nr and B — > pp 
decays. There are many high-mass states as well, and they can lead 
to additional phases (the inclusion of ttcli intermediate state makes 
5 « -40°). 



5 Taking penguins into account: shifts of Aq, 5 
and a and the values of C + _ and Coo 

Let us analyse to what changes of the parameters introduced and 
calculated in Section 2 penguins lead. Since QCD penguins contribute 
only to / = amplitude the value of extracted from B + q remains 
the same, see (ED- The requirement that the numerical values of B + 
and Bqq are not shifted when penguins are taken into account leads 
to the following shifts of the amplitude Aq and phase difference 5: 



Aq -> Aq + Aq , S -> 5 + 5 , 

v td 



5 = -V6 



Vub 
Vtd 



cos a cos(<5 p + 5q — 5q)P , 



V, 



ub 



cos a sin(5p + 5q — 5q)P/Aq 



(36) 
(37) 

(38) 



where only the terms linear in P are taken into account. For numerical 
estimates we take: 



ub 



sin 7 
sin (3 



2.3 ±0.2 



(39) 



where (3 = 22°, 7 = 60° ± 10°. In the factorization approach we have 
(see Appendix) 7 : 



Pf 



2mi 



-CI4 



-a 6 = 0.06 



[m u + m d )m b 
and shifts of Aq and 5 are small: 

- 0.12 <A < 0.12, -4°<<S<4 C 

for 



(40) 



(41) 



A = 1.5 , -1 < cos(£ p +<5 -£o) , sm(6p+5o-5o) < 1 and 70° < a < 110 c 

(42) 



7 Note that the definition of P used in the present paper differs in sign from that in 0]. 
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In particular even if the penguin contribution is underestimated by 
factor 2, the statement that S + S is large still holds 8 (note that a can 
be closer to 90°). 

The following two equations for direct CP asymmetries determine 
P and S p + Sq — So (as far as Ao, A 2 and S = S 2 — Sq are known): 



Ca 



P 



x/3 



sma[V2A sin(<5 — Sq — S p ) + A 2 sin(S 2 — Sq — S p )]/ 



Al A\ A A 2 
6 12 3x72 



cos S 



\3 



AqP cos a cos((5q — So — Sp) 



A 2 P 
V3 



cos a cos(S 2 — So — S p ) + P 2 } , 



(43) 



Coo — 



2 ~ 

* -Psina[^ sin((5o - S - S p ) - V2A 2 sin(<5 2 -6q — S p )}/ 



I [§ + § -^-A A 2 co S 6 



\ \ cos a cos(£ — 6q — S p ) + 



+ —7^A 2 P cos a cos(S 2 — So — S p ) + P 2 } , 
v3 



where 



P - 



P « 2.3P 



(44) 



(45) 



Three last terms in denominators of ( 1331) and (JHJ) lead to less than 
10% variations of the numerical values of C + _ and Coo for P < 0.3. 
Neglecting them we get: 



1 - 0.57^ 



sin(£ 2 ~ ^0 ~ £ p ) 

sin(^-^o-y ~ 1-4 + 0.41^ 



(46) 



where the numerical values for A 2 , Aq and S from (Ej), (El) and ( JT2T) 
correspondingly were used. 

From the central values in the last column in Table 1 of C + _ and 
Coo we get: 

So - So -S p = 70° , P = 0.11 . (47) 

The numerical value of P is two times larger than the factorization 
estimate of it presented in ( 1301) , while Sq — Sq — S p largely deviates from 



indication of such an underestimate follows from the probability of b — > s penguin dominated B + —> K°ir + 
decay. 
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30° which is our estimate of 5o, while Sq should be considerably smaller 
as well as 5 P the latter being close to 30° only for very asymmetric 
configurations of quarks in 7r mesons and is smaller otherwise |5j. If 
the experimental accuracy of were good we would be able to use 
the results obtained for determination of the value of the angle a from 

S.\ , realizing in this way Gronau-London approach [T]. 

However the experimental uncertainty in Coo is very big, while the 
measurements of C + _ by Belle and BABAR contradict each other. 
So let us look which values of the direct asymmetries follow from our 
formulas. 

Denominator of the expression for C H is close to one and in the 

expression in brackets in nominator first term dominates. Neglecting 
60 and 5 P and taking £0 = 30°, 62 = —20° we get C + _ = —0.04 for the 
value of penguin amplitude obtained in the factorization approach, 
(JH3J). We reproduce BABAR central value of C + if we suppose that 
factorization underestimate penguin amplitude by factor 2; however in 
order to reproduce Belle number we should accept that factorization 
is wrong by factor 10, which looks highly improbable. 

What to do if C + _ appeared to be equal to the average of the 
present day Belle and BABAR results C + _ « —0.3? One possibility 
is to suppose that 62 — So — 5p « 0, while £0 — $0 — 5p « 50° and to 
look for FSI mechanism which provides such a result 9 . 

Coo is also negative while its absolute value is larger than C + : 
denominator is about .55 while in nominator both terms are negative. 

The requirement that the value of CP asymmetry S+- is not changed 
when penguins are taken into account leads to the following shift of 
the value of the unitarity triangle angle a: 

a = a T + a , (48) 

a = i=sma[V2AQCOs(5Q — 5o — 5p) + A2Cos(52 — 5o — 5p)} . (49) 

2\/3 

Substituting the numerical values of A2 from QZj), Aq from (JHJ), P from 
( 1431) and substituting sin a by one and both cos by 0.9 we get: 

a « -2.3P « -7°, (50) 

where the result of the matrix element of the penguin operator cal- 
culation in factorization approximation (I40J) is used. We observe that 

9 Let us note that in paper |24| argument in favor of ~ —0.3 is presented which is based on the comparison 

of the direct CP violation in B(B) — > tt + tt~ and B{B) — » K + n~ (K~tt + ) decays (see also [2HJ ) - 
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our approach is at least selfconsistent: the shift of a due to penguin 
contribution is small. For the BABAR value of S + - we obtain: 

«BABAR = «BABAR + « = 92 ° ± 5 ° ■ 

In the case of the averaged experimental values we get: 

^average — ^average ^ — 98 i 4 . (<-^2) 

Theoretical uncertainty of the value of a can be estimated in the 
following way. Let us suppose that the accuracy of the factorization 
calculation of the penguin amplitude is 100% (in all the examples 
considered in this paper it was much better). Then: 

a = -7°±7° theOT , (53) 

^average = 98 =L 4 exp i 7 tlieor , (54) 

while BABAR value is smaller: 

«babar = 92° ± 5° xp ± T them . (55) 

Better theoretical accuracy of a follows from B — ► p + p~ decays, 
where penguin contribution is two times smaller. Since FSI phases are 
small in these decays, results of the paper [4J are directly applicable: 

a?"> = 92° ± 7« xp ± 4 t ° heor , (56) 

where we take the WHOLE penguine contribution as an estimate of 
the theoretical uncertainty. 

The model independent isospin analysis of B — > pp decays per- 
formed by BABAR gives |26j : 

<W = 100°±13° , (57) 

while the analogous analysis performed by Belle gives j2Zj : 

<eiie = 87° ± 17° . (58) 

Finally, the global CKM fit results are|2Bll2S|: 

acKM f mer = 97 ±5° , a UTfit = 95° ±5° . (59) 
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6 Conclusions 



1. The moduli of the amplitudes Aq and A% of B decays into tttt states 
are given with good accuracy by factorization of the tree quark 
diagram, while FSI phase shift between these two amplitudes is 
very large, \5\ « 50°, which explains large B& — > 7r°7r° decay 
probability. 

2. Theoretical uncertainty of the value of a extracted from B — > tttt 
data on S + - is at the level of few degrees. 

3. Resolution of the contradiction of Belle and BABAR experimental 
data on CP asymmetries in B — > tttt decays is very important 
both for understanding the FSI dynamics (the data on C) and 
for determination of angle a (the data on S). 

We are grateful to A.E. Bondar, A.Yu. Khodjamirian and L.B. 
Okun for useful discussions. 

A.K. was partly supported by grants CRDF RUP 2-2621-M0-04 
and RFBR 04-02-17263; M.V. was partly supported by grants RFBR 
05-02-17203 and NSh-5603.2006.2. 
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Table 1. Experimental data on B — > tttt decays. Branching ratios 

are in units of 10~ 6 . 





BABAR 


"D 11 

Beiie 


TT T71 

Heavy r iavor 
Averaging Group |7j 


B+ 


5.5 ± 0.5 


A A 1 f\ i"7 

4.4 ± 0.7 


5.0 ± 0.4 




1.17 ± 0.33 


2.3 ± 0.5 


1.45 ± 0.29 




c; o _i_ n 7 
O.o zt U. / 


O.U zt l.o 


0.0 zt U.O 


S + - 


-0.30 ±0.17 


-0.67 ±0.16 


-0.50 ±0.12 


c + 


-0.09 ±0.15 


-0.56 ±0.13 


-0.37 ±0.10 


Coo 


-0.12 ±0.56 


-0.44 ±0.56 


-0.28 ±0.39 



17 



Appendix 



In order to calculate amplitude A2 in the factorization approxima- 
tion we start from the following Hamiltonian: 

Gp 



H = -4= |K&Kdl e n [ciu<y a {l + 75)6^(1 + <y 5 )u + 

+ c 2 dj a (l + 7 5 )6w7 a (l +7 5 )w] , (Al) 
and take the matrix element of it between B lt and 7r~7r° states: 



G 



F 



V2 



\V*Kd\ e 



■I'y 



[c\ + y) < 7T |d 7a (l + 7 5 )w|0 > x 
1 



x < 7r>7 a (l + -f 5 )b\B u > +(-ci + c 2 ) < 7r>7 a (l + 75 )w|0 > x 



X < 7T |G?7 a (l + 7 5 )6|5 u > 

GF \VubV: d \ e _i7 |(ci + c 2 ) < 7T~\ J 7a (l + 75 )w|0 > x 



(A2) 



X < 7T°|w7a(l +75)&|5 U > : 



G 



F 



x/2 



K6Kk|c" fr 3(ci + P2)/^" X 



1 



X 



/ + (o)(p a s +pr) = 



Comparing the last expression with (JHD we get 

4 = 



2 , 

ci + c 2 ) = -7=(ai + ^2) = 1-35 , 



(A3) 



where ai = c\ + c 2 /3 = 1.02, a 2 = c 2 + ci/3 = 0.15 and the values 
c\ = 1.09, c 2 = —0.21 from paper [T2j were used. 

To find the amplitude Aq in the factorization approximation let us 
take the matrix element of H between and 7r + 7r~-states: 



M 



f 



and comparing it with (CD) we get 



° F : V ub V; d \e- i ~'f„mlf + (0)(c 1 + f, 



(A4) 



(A5) 



From ( IA3I ) and ( IA5I ) we obtain: 

V2 



4 = ^y(5ci - c 2 ) = 1.54 



(A6) 
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The amplitude of B d — > tt tt° decay can be calculated analogously, or 
it can be constructed from Aq and A{ with the help of 



M 



G 



F 



V2 
Gp 

72 



5ci 



V ub V: d \e-^f 7r m 2 B f + (0) [ () 
K^|e" i7 / 7 rm|/ + (0)(c 2 + ^ 



C2 8ci + 8c 2 



9 



(A7) 



a well known prediction of smallness of B d — > 7r°7r° decay probability in 
the factorization approximation. However if FSI phase shift between 
Aq and A% is large, this compensation disappears. 

QCD penguins generate an additional term to the weak interaction 
Hamiltonian, ( 1A1I ): 



AH = 



G 



F 



V2 



■)VtbV t * d [c 3 3 + C4O4 + C5O5 + c 6 6 



(A8) 



c 3 = 0.013 , 
c 4 = -0.032 , 
c 5 = 0.009 , 
c 6 = -0.037 . 



03 = dj a (l + 75)6 [uj a (l + 7 5 )w + dj a (l + 7 5 )d 

4 = d 7a(l + 75)& C [^c7a(l + 75)« a _+ d c 7a(l + 75^' 

5 = d7o(l + 75)^ [uj a (l - 75)^ + dj a (l - 75 )d] , 

6 = d 7a(l + 7s)& c [«c7a(l - l5)u a + d c j a (l - j 5 )d 
In order to find the penguin contribution to B — > 7T7t decay am- 
plitudes let us calculate the matrix element of ( IA8I) between B d and 
7r + 7r~ in the factorization approximation. For this purpose it is con- 
venient to use Fierz transformations of 7-matrices rewriting AH in 
the following form: 



A CjT r — 

AH = --j=V tb V td [a 4 w7a(l + 75)6^(1 + 7 5 )w- 

- 2a 6 w(l + 7 5 )6d(l - 75) w] , 

where (14 = C4 + 1/3C3 = —0.027, olq = cq + l/3cs = —0.034. 
Calculating the matrix element: 



(A9) 



< ir + 7r~\AH\Bd >= 



G 



F 



V2 



Vt b V td \e^f + (0)f w m 



«4 + 



2mla 6 



and comparing this expression with (JH) we get: 



[m u + m d ){m b - m u ) 
(A10) 



2mt 



-(L4 



m u + m d )m b 



a 6 = +0.06 ±0.01 



(All) 



where m u + m d = 9 ± 3 MeV, m b = 4.5 GeV were substituted. 
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